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Introduction 

The  following  study  is  an  attempt  to  better  understand  the  problem  of  scaling  water  flows  caused  by 
shallow  explosive  charges  in  field  experiments  and,  in  particular,  the  motions  of  objects  in  the  water  that  are 
swept  along  by  the  flows.  The  motivation  is  to  be  able  to  use  relatively  inexpensive  small-scale  tests  to  obtain  the 
same  information  produced  by  expensive  full-scale  tests. 

The  Navy  is  currently  interested  in  using  general-purpose  bombs  to  clear  obstacles  that  have  been 
submerged  in  shallow  water  by  adversaries  to  impede  landing  craft.  Obstacles  of  particular  interest  are  steel 
constructions  know  as  tetrahedrons  and  hedgehogs  (illustrated  in  Figure  1).  Current  research  has  determined  that 
such  obstacles  are  moved  by  an  explosion  primarily  as  a  result  of  drag  forces  created  by  the  water  flow,  although 
the  shock  wave  also  plays  a  minor  role  for  close-in  obstacles.  U.S.  efforts  to  study  the  motions  of  obstacles  using 
small-scale  explosive  tests  have  been  underway  since  1999. 


Quite  obviously,  the  Navy’s  interest  in  scaling  the  flow  of  water  has  a  long  history.  Rigorous  research  on 
scaling  the  water  flows  of  underwater  explosions  began  during  World  War  II  and  continued  for  the  next  several 
decades.  Snay1  produced  a  very  readable  and  comprehensive  account  of  this  work  that  addressed  many  topics  and 
specifically  addressed  the  scaling  of  underwater  explosion  shock  wave  and  bubble  phenomena.  Fie  showed  that, 
while  the  shock  wave  can  be  scaled  in  small-scale  field  tests,  the  scaling  of  bubble-induced  flow  phenomena  has 
always  required  special  equipment,  such  as  centrifuges  and  vacuum  tanks,  and  has  never  been  accurately 
accomplished  in  the  field.  In  this  report  we  examine  the  problem  of  scaling  bubble  effects  in  the  field  and  suggest 
a  way  by  which  this  can  be  achieved  for  the  obstacle  clearance  problem.  Our  approach  is  to  seek  out  the 
conditions  under  which  the  basic  equations  governing  the  flow  phenomena  are  invariant  under  changes  of  scale. 
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Scaling  of  Conservation  Equations  for  Incompressible  Inviscid  Flows 

We  examine  the  requirements  for  scaling  by  first  considering  the  conservation  equations  for 
incompressible  inviscid  flows.  Szymczak2  has  shown  that  explosion  bubble  phenomena  near  the  free  surface  can 
be  accurately  modeled  by  satisfying,  over  the  fluid  and  gas  domains,  the  mass  and  momentum  conservation 
equations 


A+V-(pu)  =  0 


(1) 


(pu),  +V-(puu)  =  -pgk-Vp,  (2) 

subject  to  various  constraints  and  justifiable  assumptions  regarding  the  bubble  gas,  cavitated  regions,  and  the 
ambient  atmosphere.  Other  important  contributors  are  acknowledged  in  Reference  3.  In  the  water,  the  flow  is 
treated  as  incompressible  and  the  equations  simplify  to 


Vu  =  0 


px\t  +  pV-(uu)  =  -pgk-Vp 


(3) 

(4) 


with  pressure  and  velocity  initial  conditions,  and  boundary  conditions  applied  at  the  water  surface  and 
bottom,  at  the  gas-water  surface  that  forms  the  explosion  bubble,  and  at  the  boundaries  of  cavitation 
regions,  when  they  exist.  In  axisymmetric  cylindrical  coordinates,  Equations  (3)  and  (4)  are 


1  dru  dw  .  .  .  . 

- h  —  =  0  (Continuity) 

r  dr  dz 


(5) 


P 


du  du  du 

- b  U - b  W - 

dt  dr  dz 


dp 

= - (Horizontal  Momentum) 

dr 


(6) 


P 


dw  dw  dw 

- b  it - b  w — 

dt  dr  dz 


dp 

pg  -  —  (Vertical  Momentum), 

dz 


(V) 


where  u  and  w  are  the  velocity  components  in  the  radial  and  vertical  directions. 


By  introducing  horizontal  and  vertical  characteristic  lengths  Lr  and  Lz ,  respectively,  and  a  characteristic 
time  T ,  we  can  put  the  left  sides  of  Equations  (5),  (6),  and  (7)  into  dimensionless  forms.  Upon  substitution  of  the 

Ly  ^  -^7  H4  ^ 

quantities  u  =  —u  ,  vv  =  -^-w  ,  r  =  Lrr  ,  z  =  Lzz  ,  and  t  =  Tt  ,  where  starred  terms  are  dimensionless,  we 
obtain 


*  Lr 
Lrr  — 

'  rr 


v 


T  * 

L,,r 


d\L,,r 


-w 


-  +  - 


d\Lzz 


=  0 


(8) 
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P 


di^u' 
VL  J 

din " 


L, 


'  r  * 

—  u 


di^u' 

l T  J 

d(Lrr* 


L~  : 
+  —  w 
T 


di^u' 

[t  J 

d(Lzz* 


dp 


d  I  Lrr 


(9) 


P 


Lz  ^ 
—  w 
T  j 


L 


d\  Tt 


'r  * 

—u 


U 


—  w 
v  T  J 


d  I  Lrr 


L7  * 

+  —  w 
T 


Lv 


—  w 
kT  j 


d\Lzz 


=  Pg- 


dp 


d\Lzz 


(10) 


Rearrangement  and  some  cancellations  yield  the  dimensionless  equations 


S*  SP  SP  /-N  ^ 

1  or  u  ow 

-+ —  =  o 


*  ~  * 
r  or 


dz 


/-v  *  /■>  * 

ou  *  ou 

- h  U  - 


-  +  w 


dt 


dr 


du 

a 7* 


T^]_dp 

Lr  p  dr 


OW  *  ow 
- YU  - 


-+w 


dw*  T2  ,  T2  1  dp 


dt* 


dr 


dz 


L, 


-gQg 


Lz  p  dz 


(11) 


(12) 


(13) 


It  is  noted  that  these  results  require  a  single  characteristic  time  rather  than  different  characteristic  times  for  each 
direction.  In  Equation  (13),  we  have  replaced  g  by  gQg' ,  where  we  have  chosen  the  acceleration  of  gravity  at  sea 
level,  go,  as  our  characteristic  gravity,  and  g'  is  a  dimensionless  gravity.  For  reasons  that  will  later  become  clear, 

we  write  this  as  g'  rather  than  g*  .  We  can  also  call  g'  the  gravity  scale  factor. 

The  partial  derivatives  on  the  right  hand  sides  of  Equations  (12)  and  (13)  can  be  made  dimensionless  by 

transforming  to  r  and  z  ,  and  introducing  the  characteristic  pressure  P  =  LrLz  /  pT  ,  so  that  p  =  Pp  .  In  this 
manner  we  can  express  the  conservation  equations  as 


/-s  *  *  /-s  * 

1  dr  u  dw 

- + - =  0 

*  ~  ~  * 

r  or  oz 


(14) 


OU  *  ou 


^  uii  ^  du  —  I  dp 

+  u  - +  w  - =  -N, 


dt *  dr* 


dz 


V  dr* 


(15) 


OW  *  ow 
- Y  U  - 


-  +  w 


7  ^  * 
*  dw  Z7-1  Ar  1  8P 


dt 


dr 


dz * 


■  =  F~  —  NZ 


dz * 


(16) 
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where  the  quantities  Nr ,  Nz  ,  and  F  are  given  by 

2  2 

Nr=^,  Nz=^,  and  F  =  — .  (17) 

T-P  -  T-P  g()g’T2 

Nr ,  Nz  ,  and  F  are  respectively  known  as  the  Newton  numbers  in  the  r  and  z  directions  and  the  Froude 
number.  The  Newton  numbers  are  associated  with  dynamic  inertial  effects  and  the  Froude  number  is 
associated  with  gravitational  effects. 

The  model  and  full-scale  tests  will  be  governed  by  the  same  dimensionless  equations  if  the  Newton  and 
Froude  numbers  appropriate  for  the  model-scale,  which  we  shall  denote  by  primes,  equal  those  appropriate  for  the 
full  scale.  That  is,  if  N'r  =  Nr ,  N'z  =  Nz  ,  and  F’  =  F ,  it  may  be  said  that  the  model  and  full-scale  equations 
then  possess  dynamic  similarity,  in  the  respective  directions,  and  also  gravitational  similarity.  The  scaling 
conditions  are  thus 


1 

pr 

1  T2P 

nr 

f  pF 

2  P' 

(l1  > 

V 

p' 

L'2 

P  L2 

U1 

vr) 

p  - 

v  Lr  J 

p 

1  T'2P'  1  T2P 

p'  z;2  P  L2z 


or 


f  rjir\ 


2 


P' 


p 


p 


\Lz  J 


or 


f 

V 


T) 


(18) 


(19) 


(20) 


The  time,  pressure,  length,  and  density  scale  factors  appear  on  the  right  hand  sides  of  Equations  (18),  (19) 

pf  p!  pt  1 

,  and  (20).  Defining  the  symbols  r  =  T'  IT,  (p  =  — ,  Xr  =  — ,  Xz  =  — ,  and  n  =  — ,  to  identify  these  scale 

P  Lr  Lz  P 

factors,  respectively,  we  can  rewrite  the  scaling  conditions  as 

z2(p  =  nX2  (21) 

t2  cp  =  ^>l2  (22) 

r2g'  =  4  (23) 

Equations  (21)  and  (23)  call  for  geometric  similarity,  i.e.,  Ar  =  Az  =  A  ,  which,  in  turn,  requires  equality  of  the 
Newton  numbers  Nr=  Nz  =  N . 


The  remaining  scaling  conditions  are 

(24) 

(25) 


2  -j  2 

r  (p  =  ka 

2  ,  , 

T  g  =  A. 
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Satisfying  Equations  (24)  and  (25)  in  different  ways  leads  to  different  types  of  equivalent  experimental  designs, 
all  of  which  yield  the  same  flow  information.  All  of  these  may  be  referred  to  as  Froude  and  Newton  scaling.  The 
possibilities  of  interest  are  listed  in  Table  I.  These  include  high  gravity  and  low  pressure  scaling,  which  have 
been  successfully  used  in  the  past  to  scale  explosion  bubble  phenomena,  and  also  high  density  and  the  mixed  type 
scaling  where  the  density  and  pressure  of  the  model  both  differ  from  those  of  the  full-scale  test.  High  gravity 
scaling  requires  the  use  of  a  centrifuge,  while  low  pressure  scaling  (sometimes  referred  to  as  vacuum  tank  scaling) 
requires  an  enclosed  tank  from  which  the  atmosphere  can  be  partially  evacuated.  All  scale  factor  values  in  Table 
I,  other  than  unity  and  the  pressure  scale  factory  in  the  case  of  the  mixed  design  (E5),  are  expressed  as  functions 
of  the  length  scale  factor,  which  at  this  point  is  arbitrary.  Under  E5,  (p  is  also  arbitrary. 


Table  1.  Various  Possibilities  for  Exact  Scaling  of  Incompressible  Inviscid  Flows 


Experimental 

Design 

Length 

A 

Time 

T 

Pressure 

<P 

Density 

n 

Gravity 

t 

g 

Scaling 

Description 

El 

1 

1 

1 

1 

1 

Full-Scale 

E2 

A 

A 

1 

1 

HA 

High  Gravity 

E3 

A 

\[a 

A 

1 

1 

Low  Pressure 

E4 

A 

\[a 

1 

HA 

1 

High  Density 

E5 

A 

■\[a 

(P 

<P/A 

1 

Mixed  Type 

Initial  and  Boundary  Conditions 

Equality  of  the  model  and  full-scale  dimensionless  conservation  equations  is  not  alone  sufficient  for 
achieving  flow  similarity.  In  addition,  the  dimensionless  initial  and  boundary  conditions  for  the  model  and  full- 
scale  flows  must  also  be  identical.  The  initial  condition  requires  similarity  of  the  velocity  field  \(r,z,t)  at  similar 
initial  times,  to  and  t'o  .  We,  thus,  require 


T  /  z  l()\  *  /  *  *  ,  *  ^  ^  ^  \ 

I  t’i'r  r  ,z  ,  0  =T  T’T’r  : 


(26) 


or  equivalently, 


\(r,z,t0)  =  —\'(r',z',t'0),  r  =  Ar,  z'  =  Az,  tQ  =  rtQ: 
A 


(27) 


for  all  points  in  the  fluid  domain. 

The  boundaries  of  the  model  must,  of  course,  be  geometrically  similar  to  the  boundaries  of  the  full-scale 
test.  The  numerical  conditions  that  must  be  satisfied  at  the  boundaries  are  described  by  Szymczak  in  Reference  2 
(p.5).  The  pressures  must  satisfy  Neumann  conditions  on  “wall”  boundaries  and  Dirichlet  conditions  on  nonliquid 
surfaces,  such  as  those  of  gas  bubbles  and  the  air-water  interface.  Szymczak  states  that  the  Neumann  condition, 
or  value  of  dp  /  <3n  on  the  “wall”  boundary,  implies  that  the  normal  component  of  the  flow  velocity  is  zero, 

i.e.,vn  ,  where  n  is  the  unit  normal  to  the  boundary  at  the  point  (r,  z).  The  dimensionless  condition  v*n  then 
follows  immediately.  Reference  2  also  describes  a  viscoplastic  model  suitable  for  saturated  sand  boundaries  that 
we  will  not  address  as  this  review  is  motivational  rather  than  comprehensive. 
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The  Dirichlet  conditions  require  continuity  of  the  pressure  on  the  free  surfaces.  At  the  air-water  interface 

the  dimensionless  pressure  in  the  water  p*  must  equal  the  dimensionless  atmospheric  pressure  p\  .  It  then 
follows  that 


p'a  _p'  _  *  _P  _  Pa 
P'  P'  p  P  P 

P' 

or  p'A=—pA=(ppA  .  (28) 

That  is,  the  model  atmospheric  pressure  must  equal  the  full-scale  atmospheric  pressure  multiplied  by  the  pressure 
scale  factor  at  all  times  t.  In  a  closely  related  fashion,  the  adiabatically  changing  explosion  bubble  pressure  in  the 
model  and  full-scale  tests,  p'p(t')  and  pp(t)  respectively,  must  satisfy 

PW)  =  (PPB(P )  (29) 

at  all  similar  times  t  and  t'  =  Tt ,  including  the  initial  times  to  and  t\)  .  It  also  follows  that  scaling  restricts  model 
cavitation  bubble  pressures.  If  pp  is  the  cavitation  (vapor)  pressure  in  the  full-scale  test  (assumed  constant),  the 
model  cavitation  pressure  must  be 


p'c =<p  Pe¬ 


rn 


Because  pp  is  related  to  the  detonation  pressure  of  the  explosive  and  pp  is  the  vapor  pressure  of  the 
fluid  medium,  it  is  very  difficult  to  satisfy  Equations  (29)  and  (30)  when  designs  with  a  reduced  value  of  (p  are 
used,  such  as  E3  and  E5.  As  pointed  out  by  Snay  (Reference  1,  p.50),  because  of  these  difficulties  “any  model 
test  represents  an  approximation”  and  “the  problem  of  scaling  requires  a  thorough  understanding  of  the  physics  of 
the  phenomena”  to  judge  which  effects  should  be  strictly  reproduced  in  the  model  experiment  and  which  can  be 
ignored.  The  literature  on  vacuum  tank  scaling  contains  details  of  the  techniques  used  to  control  many  of  these 
effects.4 


Application  to  Explosion  Bubbles 

When  the  conservation  equations  are  applied  to  underwater  explosion  bubbles,  the  characteristic  length 
attributed  by  Cole  to  G.I.  Taylor,  that  relates  L  to  “the  total  energy  available  after  emission  of  the  shock  wave” 
(Cole5,  p  291),  is  often  used: 


L  = 


y  pgQg ) 


1/4 


(31) 


Today,  E  is  called  the  bubble  energy.  If  we  use  Equation  (31)  to  define  an  energy  scale  factors ,  we  get 


£=E^  =  ^_(  L'^4 


P 


\L  j 


t  n  4  t 
g  =nA  g  . 


(32) 
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This  shows  that  s  is  proportional  to  A^  for  testing  design  E3  ( g'  =  1 ,  n  —  1 ),  and  to  A?  for  methods  E2 
(  g'  =  1/2 ,  n  =  1),  E4  ( g'  =  1  ,n  =  HA. ),  and  E5  (g'  =  1  ,n  =  (p/ A  ). 

In  their  review  of  spherical  bubble  theory,  Szymczak  and  Solomon1  (referencing  Friedman6 )  stated  that 
the  bubble  energy  E  is  usually  assumed  to  be  proportional  to  the  explosive  charge  weight  W  (mass).  The  constant 
of  proportionality  Q  was  shown  to  be 


Q  = 


t3 

J  cl 


(33) 


where  J ^  is  related  to  the  usual  bubble  radius  coefficient  J  by  the  equation 


J  J x> 


1  -qW-r) 
l-a~3r 


(34) 


Here  a  =  Hmax  /  Am jn  is  the  ratio  of  the  maximum  and  minimum  bubble  radii,  and  y  is  the  ratio  of  specific  heats 
of  the  bubble  gas,  which  is  assumed  to  behave  adiabatically.  Szymczak  and  Solomon  pointed  out  that  J  is  not  a 
constant  because  a  is  a  function  of  the  charge  depth  (or  hydrostatic  pressure),  but  that  may  be  regarded  as  a 
constant  characteristic  of  the  explosive.  While  this  has  been  known  (e.g.,  see  References  1  and  7),  explosives 
handbooks  and  test  reports  have  usually  tabulated  J  values  instead  of  values.  Thus,  the  energy  scale  factor 
can  be  expressed  as 


s  = 


E' 


E 


f 

u  OQ 

\A<x  j 


IT 

W 


(35) 


and  we  may  substitute  J' /  J  for  Jlr,  /  J-f  :  when  both  the  full-scale  and  model  charges  are  subjected  to  the  same 
hydrostatic  pressures.  Otherwise,  an  appropriate  correction  should  be  made  based,  for  example,  on  the  theory  of 
Szymczak  and  Solomon. 

By  combining  Equation  (35)  with  Equation  (32),  the  various  scale  factors  of  Table  I  can  be  expressed  in 
terms  of  the  charge  weights.  These  are  shown  in  Table  II.  The  scale  factors  of  Table  II  apply  to  all  flow  quantities 
of  the  model  test  and  all  of  its  boundary  and  initial  conditions.  In  the  case  of  E3,  it  should  be  noted  that,  because 
the  length  scale  factor  is  proportional  to  the  fourth  root  of  the  charge  weight  ratio,  it  cannot  be  used  to  size  the 
model  charge.  A  only  pertains  to  the  model  charge  when  the  cube  root  applies  and  when  the  same  explosive  is 
used  in  both  the  model  and  full-scale  tests. 
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Table  2.  Scale  Factors  Expressed  in  terms  of  Charge  Weights 


Length 

X 

Time 

T 

Pressure 

<P 

Density 

n 

Gravity 

f 

g 

El 

1 

1 

1 

1 

1 

E2 

V  r  1 

Jco 
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3 

W' 

w 

1/3 

1 

1 

M)V1 

J  W 

v  °o  y 

-1/3 

E3 
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Jco 

/wy 

1/4 
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3 

W' 

W 

■1/8 
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J  w 
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1/4 

1 

1 

E4 

( r  N 

Jco 

>3  1 
W' 

w 

1/3 

( r  ^ 
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3  1 
W' 

w 

1/6 

1 

J'  \ 

/®y 

i _ Z _ i 

-1/3 

1 

E5 

( r  N 

J  co 

>3  1 
W' 

w 

1/3 

(  J'  ' 

Jco 
y^CO  J 

3  1 
W' 

w 

1/6 

<P 

<P 

(  J' 

J  CO 
y^CO  , 

\3 

W' 

W 

/ 

--1/3 

1 

Predictions  of  any  full-scale  quantity,  such  as  the  horizontal  velocity  component  u  induced  by  an 
explosion  at  the  point  (r,  z)  at  time  t,  can  be  obtained  from  the  measurement  of  the  quantity  at  the  homologous 
point  (V,  z')  and  time  t'  in  the  small-scale  test  by  applying  the  appropriate  scaling  relationship.  The  scaling 
relationships  can  be  obtained  by  expressing  both  the  full-scale  and  model-scale  quantities  in  dimensionless  forms 
and  equating  them  as  was  done  for  Equation  (26).  This  follows  because  both  comply  with  the  dimensionless 
conservation  equations.  Thus,  for  the  u  velocity  component  we  have 


T 

L 


«  ( r  ,Z  ,t  )  = 


rr  t  t  ,t 

r  z  t 
—  u  ( — ,  — , — 
L'  L'  L'  T 


(36) 


Hence,  we  express  the  full-scale  quantity  u  in  terms  of  the  small-scale  measurement  u’  as 


u(r,z,t )  =  —u'(r',z',t'),  r'  =  Xr, 
X 


Z  =  X Z,  t'  =  Tt  . 


Similar  expressions  can  be  developed  for  any  other  quantities. 


(37) 


Application  to  Shallow  Small-Scale  Explosions  in  the  Field 

Small-scale  field  experiments  conducted,  for  example,  in  ponds,  lakes,  and  filled  quarries  are  problematic 
when  approached  using  traditional  bubble  scaling  methods.  Table  II  shows  that  Hopkinson  scaling,  or  the  use  of 
cube  root  length  and  time  scale  factors,  does  not  apply  to  bubble  phenomena  unless  the  gravity  scale  factor  is 
raised,  which  cannot  be  accomplished  in  the  field.  High  density  scaling  could  be  used  in  the  field  if  a  suitable 
inviscid  higher  density  fluid  could  be  found.  Low  pressure  scaling  could  be  used  in  the  field  if  a  suitable  structure 
could  be  built  over  the  body  of  water  to  allow  reduction  of  the  atmospheric  pressure.  The  Mixed  Type  scaling 
could  be  used  to  ease  the  pressure  load  on  the  structure,  but  would  again  require  the  involvement  of  special  fluids. 
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Of  the  several  methods  considered  above,  the  low  pressure  method  is  the  closest  to  satisfying  the 
requirements  for  exact  scaling  in  the  field  because  it  involves  the  same  fluid  used  in  the  full-scale  experiment  and 
is  conducted  under  normal  gravity.  In  fact,  all  of  the  requirements  for  low  pressure  scaling  can  be  met  in  the  field 
except  some  of  those  related  to  the  scaling  of  the  pressures. 

In  the  field,  the  principal  scaling  difficulty  is  caused  by  our  inability  to  control  the  static  atmospheric 
pressure.  The  hydrostatic  pressure,  due  to  the  weight  of  the  water,  is  reduced  by  the  geometric  scale  factor  and  is, 
therefore,  properly  scaled.  Hence,  the  forces  that  distort  the  bubble  shape  are  properly  scaled  when  the  water 
boundaries  and  bubble  geometry  are  properly  scaled.  Scaling  fails  in  field  experiments  primarily  because  the 
normal  atmospheric  pressure  suppresses  the  growth  of  the  bubble  and  produces  bubbles  that  are  too  small.  Indeed, 
shallow  field  explosions  that  vent  under  full-scale  conditions  often  fail  to  vent  in  small-scale  field  experiments. 
The  reduced  atmospheric  pressure,  required  by  low  pressure  scaling,  allows  the  bubble  to  grow  larger  than  it 
would  grow  under  the  normal  atmospheric  pressure.  Hence,  under  low  pressure  scaling  the  reduction  of  the 
atmospheric  pressure  is  simply  a  way  of  achieving  the  correct  bubble  size.  This  suggests  that  if  the  size  of  the 
bubble  could  be  properly  controlled,  the  rest  of  the  flow  would  behave  in  a  manner  similar  to  the  full-scale  flow. 

These  considerations  have  led  to  the  realization  that  low  pressure  scaling  rules,  with  the  exception  of  the 
pressure  rule,  might  apply  in  the  field  if  the  bubble  can  be  grown  to  the  proper  size.  The  three  methods  available 
for  growing  the  bubble  are  (1)  to  vary  the  energy  released  in  the  explosion,  (2)  to  change  the  charge  depth,  and  (3) 
some  combination  of  both.  Varying  the  charge  weight  would  change  the  amount  of  explosive  gas  in  the  bubble 
and  distort  the  bubble  collapse  and  the  flow  during  later  bubble  cycles.  Thus,  method  (1)  is  attractive  for 
modeling  bubbles  that  vent,  but  would  appear  to  be  less  attractive  for  modeling  nonventing  bubbles.  Changing 
the  model  charge  depth  to  adjust  the  bubble  size  presents  the  obvious  problem  of  changing  the  flow  geometry. 

The  consequences,  however,  may  be  insignificant  when  the  change  in  depth  is  small  relative  to  the  overall  bubble 
size.  In  particular,  method  (2)  may  also  be  effective  for  scaling  venting  explosions  in  shallow  water. 

When  venting  occurs,  the  pressure  in  the  bubble  cavity  and  the  atmospheric  pressure  rapidly  equalize. 
Indeed,  the  assumption  that  the  bubble  pressure  changes  to  the  air  pressure  instantaneously  upon  venting  is 
employed  in  the  BUB2DS  and  BUB3DS  codes  (Reference  3,  p2-4).  Thus,  any  additional  bubble  gas  associated 
with  method  (1)  above  would  simply  be  exhausted  upon  venting  and  would  not  appear  to  complicate  the 
subsequent  flow.  For  shallow  shots,  the  bubble  usually  vents  before  its  internal  pressure  reaches  the  vapor 
pressure;  hence,  it  is  also  unlikely  that  cavitation  would  be  a  complicating  factor. 

The  relationships  between  scale  factors  and  charge  weight  listed  in  Table  II  cannot  be  expected  to  be 
valid  for  flows  induced  by  venting  bubbles.  The  bubble  energy,  or  energy  available  to  drive  the  flow,  can  no 
longer  be  assumed  to  be  proportional  to  the  charge  weight.  This  is  because  differing  venting  conditions  result  in 
differing  fractions  of  the  explosion  energy  being  exhausted  to  the  atmosphere.  Instead,  the  relationships  between 
scale  factors  and  charge  weights  for  venting  bubbles  must  be  found  by  numerical  simulations  of  the  flows. 

In  collaboration  with  Szymczak  at  NRL  and  Wardlaw  at  IHDNSWC,  both  methods  (1)  and  (2)  were 
investigated  computationally  with  favorable  results.  Using  the  BUB2DS  incompressible  code,  Szymczak 
simulated  the  behavior  of  a  vertical  Mk  82  bomb  with  its  nose  touching  the  bottom  in  six  feet  of  water  for  both 
full-scale  and  l/12th-scale  field  conditions.8  In  each  case,  he  assumed  an  initial  cylindrical  bubble  with  a  volume 
equal  to  the  volume  of  the  actual  charge.  For  the  l/12th-scale  test,  he  used  a  7.3  grams  Comp  B  charge  that  was 
sized  using  the  low  pressure  scaling  conditions  with  the  appropriate  /  values.  The  length  to  diameter  ratio  of  the 
initial  bubble  for  the  model  charge  was  taken  to  be  2.3.  The  depth  of  the  model  charge  was  decreased  until  the 
model  cavity  width  versus  time  curve  matched  the  scaled  down  full-scale  curve.  Scaling  of  the  full-scale  curve 

was  done  by  multiplying  the  cavity  width  by  the  1/12  length  scale  factor  and  the  time  by  the  1  /  Vf2  time  scale 
factor,  as  required  under  low  pressure  scaling  rules.  The  center  depth  of  the  model  charge  was  2.4  inches  instead 
of  the  3.3  inches  needed  for  geometric  similarity.  Cavity  widths  were  measured  at  cells  just  above  the  rigid 
bottom.  Figure  2  shows  the  configuration  and  indicates  that  the  two  cavity  width  versus  time  curves  matched 
extremely  well. 
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Figure  2.  Full-Scale  Run  Configuration  and  Cavity  Width/Time  Curves  for  Mk  82  Bomb  (Blue) 

and  7.3-Gram  Charge  (Red) 


Figure  3.  Full-Scale  Run  Configuration  and  Cavity  Width/Time  Curves  for  Mk  84  Bomb  (Blue) 

and  28-Gram  Charge  (Red) 
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Wardlaw9  modeled  at  both  full-scale  and  l/12th-scale  the  flow  produced  by  a  vertical  Mk  84  bomb  with 
its  nose  touching  the  bottom  in  six  feet  of  water,  as  depicted  in  Figure  3.  Because  in  this  case,  approximately  one 
third  of  the  bomb  extends  out  of  the  water,  the  events  were  first  modeled  using  Wardlaw’s  Gemini  compressible 
hydrocode.*  When  the  explosive  bum  was  completed,  Wardlaw  linked  the  results  to  the  BUB2DS 
incompressible  code,  using  the  velocity  field  produced  by  Gemini,  and  continued  the  calculations  to  completion. 
The  model  charge  was  geometrically  similar  to  the  full-scale  charge  (assumed  to  be  cylindrical)  in  the  vertical 
direction,  but  its  diameter  was  reduced  to  achieve  the  desired  explosive  weight.  This  method  was  chosen  to  keep 
the  same  fractions  of  the  explosive  charge  above  and  below  the  water  surface  for  all  runs.  The  charge  weight  was 
varied  from  that  required  by  low  pressure  scaling  (36  gms)  until  the  cavity  width  versus  time  curve  of  the  model 
matched  that  of  the  scaled  down  full-scale  calculation.  The  best  match,  which  is  shown  in  Figure  3,  occurred  for 
a  charge  weight  of  28  grams.  That  the  quality  of  the  match  is  somewhat  below  that  shown  in  Figure  2  is  thought 
to  be  due  to  the  fact  that  the  calculation  was  more  complex.  To  obtain  an  accurate  simulation  of  the  explosive 
bum  and  to  get  the  shock  wave  well  beyond  the  bubble,  the  Gemini  run  had  to  be  rezoned  three  times.  The 
difficulty  of  controlling  the  rezoning  times  may  have  introduced 

differences.  There  was  also  a  discontinuity  in  the  linkage  of  the  codes  and  a  non-physical  pressure  pulse  caused 
by  the  incompressibility  of  the  BUB2DS  code.9 

The  finding  that  the  charge  size  needed  to  grow  the  bubble  was  less  than  that  predicted  by  low  pressure 
scaling,  rather  than  more,  deserves  comment.  It  is  thought  that  this  is  because  the  diameter  of  the  model  charge 
was  smaller  than  that  of  the  scaled  down  full-scale  charge.  Because  the  charge  length  was  kept  geometrically 
similar,  a  diameter  smaller  than  the  geometrically  similar  diameter  had  to  be  used  to  achieve  charge  weights  on 
the  order  of  those  required  by  l/4th  root  low  pressure  scaling.  The  smaller  charge  diameter  and  the  higher 
pressure  at  the  surface  resulted  in  a  more  efficient  coupling  of  the  energy  of  the  model  charge  to  the  water. 
Because  the  full-scale  charge  diameter  was  larger,  the  explosion  gases  could  more  easily  escape  the  cavity.  The 
complete  explanation  is  complicated,  however,  because  the  stronger  shock  wave  of  the  full-scale  charge  relative 
to  the  model  charge  (shock  wave  phenomena  follow  cube  root  scaling  rules)  would  probably  result  in  more 
spalling  at  the  surface.  It  is  presently  unclear  whether  this  would  enhance  or  suppress  the  growth  of  the  full-scale 
cavity. 


Figure  4  shows  selected  frames  from  movies  of  the  computations  for  the  Mk  84  bomb  and  the  28-gram 
charge  at  approximately  similar  times.  The  flows  appear  to  be  quite  similar,  although  irregularities  are  evident 
that  may  be  due  to  the  complications  mentioned  by  Wardlaw  and  listed  above.  Radial  and  vertical  components  of 
the  flow  velocity  were  recorded  at  mid-depth  stations  at  several  distances  from  the  origin.  These  are  shown  in 
Figure  5.  The  horizontal  velocity  components  appear  be  a  somewhat  better  match  than  the  vertical  components. 
Nevertheless,  the  overall  agreement  between  the  model  and  full-scale  flows  is  again  quite  good.  Moreover,  the 
agreement  between  the  velocity  components  appears  to  be  best  at  early  times. 


*  Compressible  flow  equations  follow  Hopkinson  (high  gravity)  scaling.  No  attempt  was  made  here  to  scale  the 
compressible  flow.  Rather,  the  hydrocode  was  simply  being  used  to  establish  initial  conditions  for  the  incompressible  flow. 
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Figure  4.  Mk  84  Bomb  (Left)  and  28-Gram  Charge  (Right)  Development  Sequences  at 

Approximately  Similar  Times 
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Figure  5.  Radial  and  Vertical  Components  of  Flow  Velocity  for  Mk  84  Bomb 

(Blue)  and  28-Gram  Charge  (Red) 

The  calculations  suggest  that  the  entire  velocity  field  is  matched  when  the  cavity  width  versus  time  curves 
are  in  agreement,  particularly  after  venting  time.  When  the  bubbles  have  vented  and  the  entire  free  surface, 
including  the  cavity  wall,  is  at  atmospheric  pressure,  the  velocity  field  alone  determines  the  future  history  of  the 
flow.  The  fact  that  the  model  is  at  atmospheric  pressure,  rather  than  the  reduced  level  prescribed  by  the  low 
pressure  scaling  rules,  does  not  affect  the  flow  after  venting  because  the  conservation  equations  depend  on  the 
pressure  through  the  gradient  only.  Hence,  it  appears  that  the  agreement  of  the  velocity  fields  following  venting 
guarantees  the  similarity  of  the  flows. 
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Figures  6  and  7  show  that  it  is  not  always  possible  to  match  the  cavity  width  history  by  increasing  the 
model  charge  weight  at  the  geometrically  scaled  full-scale  charge  depth.  In  this  numerical  experiment  performed 
by  Wardlaw  using  the  BUB2DS  code  (initialized  by  Gemini),  conditions  were  similar  to  that  shown  in  Figure  3 
except  that  the  water  depth  was  increased  to  12  feet.10  The  Mk  84  bomb  was  then  fully  submerged.  Cavity  width 
plots  for  48  and  42  grams  of  Comp  B,  shown  in  Figure  6,  suggested  that  a  smaller  charge  weight  was  required  to 
match  the  scaled  down  full-scale  curve.  The  charge  weight  required  by  the  low  pressure  scaling  law  was  36 
grams.  The  run  at  39  grams,  however,  showed  both  (1)  that  a  smaller  size  charge  was  needed  to  be  on  the  scaled 
full-scale  curve  and  (2)  that  the  bubbles  of  such  smaller  charges  would  not  vent.  Figure  7  shows  the  39-gm 
bubble  developing  downward  and  upward  jets  at  a  scaled  time  after  the  full-scale  bubble  had  vented.  Apparently, 
the  high  atmospheric  pressure  in  the  model  test  (relative  to  that  of  the  low  pressure  scaling  rule)  prevented  the 
bubble  from  opening. 
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Figure  6.  Cavity  Widths  Comparison  -  All  Charges  at  Depth  Geometrically 

Similar  to  Full-Scale  Depth 


Figure  7.  Frames  Showing  1/1 2-Scale  Bubble  Developing  Upward  and 
Downward  Jets  After  Full-Scale  Bubble  Has  Vented 
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Figure  8  shows  that  it  was  possible  to  match  the  cavity  width  versus  time  curve  for  the  Mk  84  bomb  in  12 
feet  of  water  by  reducing  the  depth  of  the  model  charge.10  The  charge  depth  was  changed  from  7  inches  to  4.2 
inches  to  obtain  a  charge  depth  to  maximum  bubble  radius  ratio  of  0.2  for  the  36-gram  charge,  which  was  the 
value  of  the  ratio  in  Szymczak’s  simulation  of  the  Mk  82  bomb  in  6  feet  of  water.  A  tested  empirical  rule  is  that 
bubbles  shallower  than  a  quarter  of  a  maximum  bubble  radius  vent  under  normal  atmospheric  pressure;  hence, 
this  depth  was  chosen  to  ensure  venting.  The  size  of  the  maximum  bubble  radius  was  obtained  from  the  empirical 
bubble  equation.  The  39-gram  charge  at  depth  of  7  inches  had  a  ratio  of  about  one  third,  and,  as  seen  above,  did 
not  vent. 

From  Figure  8  we  see  that  the  cavity  width  versus  time  curve  of  the  32-gram  charge  was  in  reasonably  good 
agreement  with  the  scaled  full-scale  curve.  The  36-gram  curve  is  also  shown.  Figure  9  shows  the  radial  and 
vertical  velocity  components  at  selected  horizontal  stations  for  all  three  runs.  The  radial  components  are,  perhaps, 
in  better  overall  agreement  than  the  vertical  components.  Fine-scale  differences  between  all  three  computations 
are  apparent  in  the  movies  of  the  flows,  and  these  may  be  responsible  for  some  of  the  differences  shown  in  Figure 
9.  On  the  whole,  the  agreement  between  the  model  and  full-scale  velocity  components  is  quite  acceptable. 


Figure  8.  Cavity  Widths  Comparison  -  Model  Charges  at  Depth  of  4.2  Inches  to  Ensure  Venting 
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Figure  9.  Radial  and  Vertical  Components  of  Flow  Velocity  for  Mk  84  Bomb  (Blue)  in  12-ft  of 
Water  and  32-gm  (Black)  and  36-gm  (Red)  Comp  B  Charges  at  4.2  in. 
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Scaling  of  Obstacle  Displacements 

We  now  examine  whether  obstacles  placed  in  the  small  scale  flow  field  behave  in  manners  that  are 
similar  to  their  behavior  under  full-scale  conditions.  We  examine  first  the  case  of  an  obstacle  that  translates 
without  rotation.  Computer  generated  animations  of  obstacle  motions  made  by  Wardlaw  (using  a  rigid  body 
model  of  the  motions  of  a  tetrahedron  coupled  to  Szymczak’s  BUB2DS  incompressible  hydrocode  that  has  been 
started  using,  as  in  the  case  of  the  Mk  84  bomb,  the  GEMINI  hydrocode)  showed  that,  in  some  cases,  obstacles 
did  not  rotate  appreciably  while  being  transported  by  the  flow. 1 1  Generally,  however,  obstacles  do  rotate  during 
translation.  We  will  comment  at  the  conclusion  of  this  section  on  the  effects  rotations  may  have  on  the  motion  of 
the  obstacle  and  the  issue  of  scaling. 

If  the  translating  obstacle  is  lifted  and  translated  by  the  flow,  a  first  approximation  of  the  equations 
governing  the  obstacle’s  motion  can  be  assumed  to  be 


in  the  horizontal  direction  and 


1  2 

(m  +  m )r  =  ~P  Crar  (u  -  r) 


(38) 


(m  +  mA  )z 


2 


pCzaz(w- 


mg()g' 


(39) 


in  the  vertical  direction.  Here,  m  is  the  mass  of  the  obstacle,  m Ar  and  m^z  are  added  masses  for  flows  in 
the  r  and  z  directions,  respectively,  ar  and  az  are  the  areas  of  the  obstacle  presented  to  the  horizontal  and 
vertical  components  of  the  flow,  andCr  andCz  are  the  drag  coefficients.  Equations  (38)  and  (39)  are 
based,  in  part,  on  an  analysis  performed  by  Goeller  and  Ruben,12  who  studied  the  case  of  an  obstacle 
that  is  driven  horizontally  by  the  water  flow  and  is  resisted  by  frictional  contact  with  the  bottom.  They 
concluded  that  a  “constant  drag  coefficient  appears  reasonable;  that  is,  there  is  little  dependence  on 
Reynolds  Number.”  Hence,  we  shall  assume  that  Cr  and  Cz  are  unvarying  constants. 

The  need  for  added  masses  in  the  governing  equations  is  a  currently  unresolved  issue.  They  are  not 
employed  by  Wardlaw  in  his  hydrocode-coupled  rigid  body  code  because  the  obstacle  mass  is  being  driven  by  the 
flow  rather  than  driving  it.11  We  shall  carry  them  along,  nevertheless,  in  the  event  that  they  prove  to  be  necessary 
during  some  phases  of  the  motion. 

The  obstacle  masses  and  the  area  presented  to  the  horizontal  flow  component  will  vary  if  the  length  scale 
L  of  the  obstacle  is  changed.  We  can  assume  that  the  variation  of  the  mass  of  the  obstacle  with  L  can  be 
represented  as 


m  Pobs  ^  fm  •  (40) 

Here  pQ^s  is  the  density  of  obstacle  material  and  fm  =  V //f ,  is  a  dimensionless  form  factor  that  is  independent 

of  scale,  where  V  is  the  obstacle  volume.  Similarly,  the  variation  with  L  of  the  area  presented  to  the  horizontal 
flow  can  be  written  as 


=  LZfa 


(41) 
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where  fa  is  a  dimensionless  horizontal  area  form  factor.  In  like  manner,  we  write  the  presented  area  in 
the  vertical  direction  as 


«z=£24-  (42) 

The  added  masses  of  the  obstacle  are  somewhat  more  involved.  Assuming  potential  flow,  Pozrikidis  has 
shown13  that  the  kinetic  energy  K  of  the  fluid  for  a  general  motion  of  a  rigid  body  of  volume  V  can  be  written  as 

K  =  ^pVAijWiWj,  (43) 

where  W  =  {u,v,w,  cox ,  (Ov ,  a>z  )'  is  a  vector  of  translational  and  rotational  velocity  components  of  the  rigid  body 

an&Ajj  is  the  symmetric  six-by-six  grand  added  mass  matrix.  Pozrikidis  states  that  the  “value  of  A  depends 

exclusively  upon  the  instantaneous  body  shape  and  orientation  and  presence  of  other  objects,  but  is  independent 
of  the  body’s  linear  or  angular  velocity  or  acceleration.”  The  important  idea  for  us  is  that  A  depends  on  body 
shape,  but  not  its  scale.  Since  we  are  assuming  translational  motion  in  two  directions  only,  A  has  only  three 
nonzero  terms  corresponding  to  the  r  and  z  directions  and  a  cross  term.  Because  the  volume  dimensions  are 
proportional  to  their  associated  scale  factors,  Equation  (43)  shows  that  it  is  possible  to  represent  the  virtual  masses 


by  the  forms 

mAr  =  P1-2 fmAr 

(44) 

and 

mA7  ~  PL  fmA_  j 

(45) 

where  fm  and  fm  are  dimensionless  form  factors  that  involve  elements  of  the  grand  added  mass  matrix. 

Ar  Az 

If  we  insert  Equations  (40),  (41),  (42),  (44),  and  (45)  into  Equations  (38)  and  (39),  along  with  the 

expressions  r  =  -r-r* ,  z  =  -^-z*  ,  r  =  —  r* ,  z  =  —  z*  and,  as  before, 

T2  T2  T  T 

~J~ j  ~J~ j 

u  =  —u  and  w  =  —  w  ,  in  which  the  starred  quantities  are  dimensionless,  we  get 


(Pobsfm  +  PLir  A  r*  =  \p fa,. <f  ~ f  “ >2 


(46) 


(, Pobsfm  +PfmAi  >P  4"  A-pCfi 4  (3>*  Ai*)2  ~pobsL3fmgog'  (47) 


If  we  insert  — —  =  F  gQg'  from  Equations  (17)  and  rearrange,  we  obtain 

rji  A 
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and 


..*  9  Crpfa 

r  =  z 


( Pobsfm  Pfm  A  ) 


*  *9 

(u  -r  )2 


z*  =  — 


9  CZF pfa  (w  2  )  Pobsfm 


Pobsfm  Pfm ) 


(48) 


(49) 


Equations  (48)  and  (49)  show  that  the  displacement  equations  are  invariant  under  Froude  scaling,  i.e.,  when  the 
Froude  number  F  is  the  same  for  both  model  and  full-scale  tests. 


Obstacle  rotations  will  add  one  to  three  additional  equations  of  motion,  depending  upon  the  initial 
obstacle  symmetries,  or  lack  of  symmetries,  relative  to  the  directions  of  water  flow.  Fack  of  symmetry  may  also 
create  transverse  translational  motion  and  a  three  dimensional  trajectory.  It  seems  reasonable  to  expect  that  the 
equations  of  rotation  and  translation  would  be  largely  uncoupled,  at  least  to  a  first  approximation.  Without  going 
into  a  detailed  analysis,  the  author  suspects  that  the  primary  effect  of  obstacle  rotations  would  be  to  cause  the 
form  factors  appearing  in  Equations  (48)  and  (49)  to  vary  with  time.  Since  these  are  scale  independent,  the 
motion  of  the  rotating  and  translating  obstacle  will  likely  also  be  invariant  under  scale  changes.  This  conjecture 
could  best  be  tested  with  the  use  of  Wardlaw’s  hydrocode-coupled  rigid  body  code,  which  allows  six  degrees  of 
freedom. 


Experimental  Tests  of  the  Shallow  Explosion  Scaling  Theory 

In  the  summer  of  2005,  field  tests  were  conducted  by  SRI  International,  under  the  direction  of  P.Gefken, 
using  1/12-scale  tetrahedron  obstacles  and  a  water  depth  of  six  inches.14  Comp  B  charges  weighing  7.3  grams  and 
28  grams  were  prepared  that  modeled  the  Mk  82  and  Mk  84  bombs,  respectively,  and  were  configured  in  a 
manner  consistent  with  the  theory  and  calculations  discussed  above.  Six  replicated  shots  were  conducted  for  both 
the  Mk  82  scaled  conditions  and  the  Mk  84  scaled  conditions.  For  both  cases,  all  tetrahedrons  in  three  of  the 
shots  were  placed  in  a  “toward”  orientation  and  in  an  “away”  orientation  in  the  other  three  shots.  A  toward 
orientation  means  that  the  ground  triangle  (the  side  of  the  tetrahedron  in  contact  with  the  ground)  pointed 
symmetrically  in  the  direction  of  the  charge.  The  ground  triangle  pointed  away  from  the  charge  for  the  away 
oriented  tetrahedrons.  The  displacements  of  54  tetrahedron  obstacles  were  measured  for  each  case. 

Figure  10  shows  the  radial  displacements  of  the  obstacles  in  the  7.3-gram  tests  as  a  function  of  their  initial 
radial  standoffs  (measured  to  the  point  on  the  ground  below  the  top  vertex).  Diamond  symbols  denote  the  toward 
orientations  and  square  symbols  denote  the  away  orientations.  Different  colors  are  used  to  identify  different 
shots.  The  black  solid  line  is  the  mean  obstacle  displacement  from  a  regression  fit  to  the  scaled  down  full-scale 
Mk  82  test  data,  which  are  displayed  here  as  small  red  plus  symbols.  The  two  dashed  lines  represent  ±2cr  curves 
for  the  full-scale  displacements.  In  the  analysis  of  the  full-scale  data,  the  standard  deviations  were  assumed  to  be 
proportional  to  the  mean.  This  assumption  is  not  fully  realistic,  but  it  is  in  reasonable  agreement  with  the  full- 
scale  data.  According  to  theory,  about  98  percent  of  the  displacement  points  should  fall  between  the  dashed  lines. 
The  full-scale  analysis  was  based  on  the  displacements  of  a  total  of  16  tetrahedron  obstacles  in  four  Mk  82  bomb 
tests. 


Figure  10  suggests  that  the  SRI  radial  displacement  data  agree  reasonably  well  with  the  scaled  down  full- 
scale  displacements.  At  7  inches,  the  SRI  data  lie  somewhat  above  the  full-scale  mean,  while  at  12  and  15  inches 
they  appear  to  fall  below  the  mean  curve.  The  continuing  upward  trend  shown  by  the  SRI  displacements  at  close- 
in  initial  standoffs  is  probably  realistic,  since  the  downturn  of  the  full-scale  mean  is  based  upon  only  two  data 
points.  In  theory,  one  would  expect  a  downturn  at  some  point  as  the  initial  standoff  is  reduced  to  zero.  However, 
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in  reality,  the  obstacle  may  break  up  before  this  is  observed.  The  small-scale  displacements  appear  to  fall  off 
more  sharply  with  initial  standoff  than  the  full-scale  data.  The  scatter  of  the  SRI  data  looks  much  like  the  full- 
scale  scatter,  although  there  may  be  a  local  increase  around  15  inches.  There  is  no  detectable  difference  between 
the  displacements  of  obstacles  with  toward  and  away  orientations. 

Clearly  the  assumption  that  the  standard  deviation  is  proportional  to  the  mean  is  incorrect  since  there  is 
appreciable  spread  in  the  data  after  the  mean  has  reached  zero.  A  more  correct  assumption  may  be  that  the 
standard  deviation  is  proportional  to  the  mean  of  the  total  distance  traveled.  Calculations  made  by  Wardlaw  with 
a  rigid  body  code  for  obstacle  motions  coupled  to  his  Gemini  hydrocode  show  that  the  obstacles  move  back  and 
forth  with  the  flow  to  a  considerable  extent  before  coming  to  rest.11  This  apparently  produces  the  scatter  at  the 
larger  distances.  Although  direct  observation  of  the  total  travel  distance  would  probably  be  difficult  for  each 
obstacle  tested,  suitable  information  could  probably  be  obtained  from  Wardlaw’s  codes  to  revise  the  statistical 
analysis. 

The  angular  displacements  for  the  7.3-gram  shots,  shown  in  Figure  11,  suggest  increased  scatter  for 
obstacles  placed  at  the  18-  and  21 -inch  standoffs.  This  could  also  reflect  a  dependence  on  the  total  travel 
distance.  There  does  not  appear  to  be  an  effect  due  to  obstacle  orientations,  but  more  tests  would  be  required  to 
detect  such  differences. 

A  possibly  significant  difference  between  the  full-scale  and  l/12th-scale  conditions  is  that  the  masses  of 
the  tetrahedrons  in  SRI’s  model  tests  appear  to  have  been  10  to  20  percent  low.  The  weight  calculated  for  an 
idealized  full-scale  tetrahedron  measuring  56  inches  from 
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INITIAL  RADIUS  (INCHES) 


Figure  10.  Comparison  of  7.3-Gram  Test  Radial  Displacements  with  Scaled  Full-  Scale  Mean 
(Solid)  and  ±2cr  Standard  Deviation  (Dashed)  Curves.  Colors  Denote  Different  Shots.  Diamonds 

and  Squares  Denote  Toward  and  Away  Orientations. 
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Figure  11.  Angular  Displacements  in  7.3-Gram  Tests.  Colors  Denote  Different  Shots. 
Diamonds  and  Squares  Denote  Toward  and  Away  Orientations. 
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vertex  to  vertex,  with  a  5/8-inch  thickness,  leg  widths  of  4  inches  in  both  directions,  and  using  a  steel  density  of 
0.284  pounds  per  cubic  inch,  is  418  pounds.  The  density  used  corresponds  to  a  linear  weight  of  15.7  lbs/ft  for  the 
angle  steel.  This  result  is  somewhat  lower  than  that  estimated  by  Goeller  and  Rubin  in  Reference  13.  At  1/1 2th- 
scale,  the  calculated  tetrahedron  mass  would  be  1 10  grams.  The  masses  of  the  tetrahedrons  used  in  the  SRI  tests 
were  93  ±1  grams.14  It  is  unclear  how  the  lighter  mass  might  have  affected  the  results  of  the  tests.  Previous 
experiments  conducted  by  SRI  using  50-gram  Comp  B  charges  in  water  six  feet  deep  and  involving  like  obstacles 
of  different  masses  showed  that  the  heavier  ones  tended  to  become  airborne  more  often  than  lighter  obstacles.15  It 
was  thought  by  some  participants  in  the  effort  that  the  additional  inertia  made  it  easier  for  the  heavier  obstacles  to 
exit  the  water.  It  is  reasonable  to  suspect  that  the  lower  inertias  of  the  tetrahedrons  may  have  resulted  in 
premature  stopping  and  perhaps  a  steepening  of  the  slope  of  the  mean  curve. 

Just  prior  to  publication  of  this  report  it  was  learned16  that  the  Mk  82  bombs  used  in  the  four  full-scale  tests 
had  explosive  weights  that  were  about  9  percent  below  the  standard  weight  of  192  lbs  upon  which  the  scaled  tests 
were  designed.  It  was  also  learned  that  the  explosive  used  in  the  full-scale  tests  was  H-6  instead  of  Tritonal, 
which  had  been  assumed.  Because  the  /  value  for  H-6  is  slightly  below  that  of  Tritonal  (about  2  percent  less),  the 
full-scale  displacements  in  Figure  10  represent  a  bomb  with  the  same  explosive  and  a  weight  that  is,  perhaps,  10 
percent  below  the  standard  weight.  This  difference  may  account  for  the  close-in  disagreement  between  the  model 
and  full-scale  displacements.  It  is  reasonable  to  expect  that  displacement  curve  for  the  standard  weight  would  be 
proportionately  raised  and  the  difference  between  the  slopes  of  the  full-scale  and  model-scale  data  lessened.  The 
implications  of  this  charge  weight  discrepancy  will  be  further  investigated  in  2006. 

The  results  for  the  28-gram  tests  conducted  by  SRI  are  shown  in  Figures  12  and  13.  No  full-scale  tests  of 
the  Mk  84  bomb  in  water  six  feet  deep  with  obstacles  have  been  conducted;  hence,  no  full-scale  comparisons  can 
be  made.  The  28-gram  test  results  again  show,  in  both  Figures  12  and  13,  an  increase  in  the  scatter  for  the 
obstacles  at  intermediate  distances  from  the  charge.  Qualitatively,  the  results  are  much  like  those  of  the  7.3-gram 
shots.  An  interesting  detail  seen  in  Figure  12,  although  not  in  Figure  10,  is  that  the  displacements  of  obstacles 
with  toward  orientations  (diamonds)  appear  to  be  slightly  higher  than  those  with  away  orientations  (squares)  at 
close  in  positions.  This  might  be  expected  from  a  theoretical  viewpoint  because  the  toward  orientation,  with  two 
legs  tilted  toward  the  explosion,  should  generate  more  lift  in  the  outflow  direction  and  less  in  the  inflow  direction, 
thus  favoring  a  greater  net  displacement. 
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Figure  12.  Radial  Displacements  in  28-Gram  Tests.  Colors  Denote  Different  Shots.  Diamonds 

and  Squares  Denote  Toward  and  Away  Orientations. 
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Figure  13.  Angular  Displacements  in  28-Gram  Tests.  Colors  Denote  Different  Shots.  Diamonds 

and  Squares  Denote  Toward  and  Away  Orientations. 
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Conclusions  and  Recommendations 

1.  It  appears  that  obstacle  motions  produced  by  full-scale  explosions  can  be  accurately  reproduced  in  small- 
scale  experiments  when  venting  of  the  explosion  bubble  occurs.  Model  charge  conditions  that  bring 
about  a  matching  of  the  model-scale  and  scaled  down  full-scale  cavity  width  versus  time  curves  appear  to 
successfully  reproduce  all  hydrodynamic  details  of  the  water  flow  as  well  as  the  displacements  of  scaled 
obstacles  placed  in  the  flow.  This  matching  process  can  be  achieved  by  using  hydrocodes  that  model  the 
explosion  flow  effects.  The  appropriate  time  scale  factor  is  the  square  root  of  the  length  scale  factor, 
which  can  be  set  arbitrarily  for  the  model  experiment.  Conditions  for  this  computer-guided  Froude 
scaling  are  close  to  those  used  for  low  pressure  scaling,  but  they  do  not  require  reduction  of  the 
atmospheric  pressure.  Hence,  they  are  appropriate  for  tests  conducted  in  the  field. 

2.  While  the  discussion  in  this  report  has  been  focused  on  scaling  the  flow  effects  of  a  single  explosive 
charge,  it  would  appear  that  the  method  can  be  extended  to  the  multiple-charge  scenario.  It  is  likely  that 
sizing  the  various  charges  could  be  done  on  a  one-at-a-time  basis  using  the  method  summarized  above, 
but  this  suggestion  should  be  the  subject  of  future  research. 

3.  It  is  recommended  that  the  relationship  between  obstacle  displacement  and  obstacle  mass  be  investigated 
computationally  using  Wardlaw’s  Gemini-coupled  rigid  body  obstacle  displacement  code. 

4.  It  is  recommended  that  additional  testing  be  performed  by  SRI  International  using  obstacles  with  masses 
that  closely  match  those  of  the  scaled  down  full-scale  obstacles  for  the  purpose  of  validating  the 
computer-guided  scaling  approach  discussed  in  the  report.  It  is  recommended  that  the  experimental  study 
also  explore  other  factors  that  might  affect  the  shape  of  the  mean  displacement  versus  initial  standoff 
curve,  such  as  sand  characteristics. 

5.  It  is  recommended  that  the  scaling  bubbles  that  vent  during  the  second  cycle  of  expansion  be  explored 
computationally. 

6.  It  is  recommended  that  theoretical  studies  be  funded  to  rigorously  establish  the  concept  of  computer- 
guided  scaling  of  shallow  explosions. 

7.  It  is  recommended  that  the  comparison  between  the  full-scale  and  model  tests  be  reviewed  and 
consideration  given  to  additional  small-scale  tests  designed  using  the  actual  full-scale  test  charge  weights. 
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